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SUPERSONIC FLOW PAST OSCILLATING AIRFOILS INCLUDING 
NONLINEAR THICKNESS EFFECTS 1 

By Melton D. Van Dyke 


SUMMARY 

A solution to second order in thickness is derived for har- 
monically oscillating two-dimensional airfoils in supersonic 
flow. For slow oscillations of an arbitrary profile , the result is 
found as a series including the third power of frequency. For 
arbitrary frequencies, the method of solution for any specific 
profile is indicated, and the explicit solution derived for a single 
wedge. 

Nonlinear thickness effects are found generally to reduce the 
torsional damping, and so to enlarge the range of Mach numbers 
within which torsional instability is possible. This destabiliz- 
ing effect varies only slightly with frequency in the range in- 
volved in dynamic stability analysis, but may reverse to a 
stabilizing effect at high -flutter frequencies. Comparison with 
a previous solution exact in thickness suggests that nonlinear 
effects of higher than second order are practically negligible. 

The analysis utilizes a smoothing technique that replaces 
the actual problem by one involving no kinked streamlines. 
This stratagem eliminates all consideration of shock waves 
■from the analysis, yet yields the correct solution for problems 
that actually contain shock waves. 

INTRODUCTION 

As linearized supersonic-flow theory is increasingly ap- 
plied to problems of unsteady motion of lifting wings, the 
results are sometimes advanced with the warning that they 
may be significantly affected by nonlinear effects of thickness. 
Such caution is justified because it is known that even for 
steady flow linearized theory is often inadequate for predict- 
ing the pitching moment — and prediction of moments is one 
of the main objectives of unsteady-flow theory. It may be 
anticipated that nonlinear effects will become increasingly 
important as the Mach number falls toward unity, particu- 
larly for slow oscillations. 

In the present work the effects of thickness are determined 
for a harmonically oscillating two-dimensional airfoil by 
calculating the second-order solution. This is the counter- 
part for unsteady motion of the well-known steady-flow 
result of Busemann (ref. 1). First, for slow oscillations a 
solution is found for an airfoil of arbitrary profile. The 
result is given as a series that includes terms up to the third 
power of the frequency. Second, for arbitrarily high fre- 
quencies it is shown that a solution can be found for any 
specific airfoil, and the solution is carried out explicitly for a 


single wedge. Finally, comparison is made with a previous 
solution for the wedge that is exact with respect to thickness 
(refs. 2 and 3), in order to assess the effects of nonlinear 
terms of higher than second order. 

Extensive use is made of a smoothing technique, which 
replaces the actual problem by one having no kinked stream- 
lines. This stratagem, which has been used previously and 
may prove useful in future problems, e limin ates all considera- 
tion of shock waves from the analysis. It, nevertheless, 
leads to the correct second-order solution for the actual 
problem, which does involve shock waves. 

METHOD OF ANALYSIS 

STATEMENT OF PEOBLEM 

Consider a sharp-nosed airfoil flying through still air at a 
uniform supersonic velocity and executing small harmonic 
oscillations. We shall be concerned with calculating the 
instantaneous pressure at the surface and, hence, the un- 
steady lift and pitching moment. If oscillations in the 
flight direction are neglected, a rigid airfoil possesses two 
degrees of freedom. The oscillation can therefore be re- 
garded as compounded of a rotation (pitching) and a vertical 
translation (plunging), which are not generally in phase. 

Although the iteration procedure to be employed yields a 
formal result for any Mach number greater than unity, the 
solution probably breaks down when the flow becomes sonic 
at any point. Since this occurs at a Mach number somewhat 
higher than that for bow-wave detachment, the upper and 
lower surfaces of the airfoil operate independently in the 
probable range of validity of the solution. It is. therefore 
sufficient to consider only the half field of flow lying above 
the airfoil, and this viewpoint will be adopted henceforth. 

It is convenient to seek a solution to second order in the 
airfoil thickness, but to only first order in the amplitude of 
oscillation. This is sufficient because second-order terms 
in the oscillation, although affecting local pressures, have 
no effect upon lift or moment since they are equal on the 
upper and lower surfaces. Then the pitching and plunging 
components of the oscillation can be treated separately, and 
the results superimposed. Furthermore, it is enough to 
consider only pitching about an arbitrary pivot, because the 
plunging case can be recovered by letting the pivot recede to 
infinity and the pitching amplitude tend to zero, their 
product remaining finite. Thus, from the point of view of an 


i Supersedes NAOA TN 2882, "Supersonic Flow Past Oscillating Airfoils Including Nonlinear Thickness Effects” by Milton D. Van Dyke, 1953. 
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observer moving with the mean speed of the airfoil (or testing 
it in a wind tunnel) the airfoil is exposed to a uniform super- 
sonic stream and oscillating slightly about a fixed pivot 

(fig- 1)- 


y 

\ 

■ 



Fiqure 1. — Coordinate system for oscillating airfoil. 


Choose the origin of coordinates at the mean position of 
the leading edge, with the x axis extending in the direction 
of the free stream. Then it is convenient to describe the 
upper surface of the airfoil in its mean (zero angle of pitch) 
position by 

y=7(x)=tf(x) (1) 

All symbols are defined in Appendix A. Here 6 is a small 
parameter representative of the airfoil thickness, so that the 
function / is of order unity. Now let the airfoil pivot about 
a point on the x axis lying a distance b downstream from the 
leading edge, and perform harmonic oscillations of frequency 
w and amplitude 0 O , so that the angle of pitch,* which is the 
angle between the instantaneous position of the airfoil and 
its original mean position (fig. 1) is given by 

0=0 O cos Gjf=0 o e 1 "‘ (2) 

(Here, as in all that follows, it is implied that actual physical 
quantities are given by the real parts of their complex repre- 
sentations.) Then, at any instant the moving upper surface 
of the airfoil is described by 

y=ef(x)—0 o e tat (x—b) (3) 

with an error of order (6*0, e0*), which is of third order and, 
consequently, negligible in the present second-order analysis. 

PERTURBATION EQUATION 

The entropy changes due to shock waves are of third order 
in the airfoil thickness and angle of pitch. Hence, to second 
order the flow is irrotational and isentropic. Because it is 
irrotational, there exists a potential function G whose 
gradient yields the velocity vector: 

-♦ u=Q I ) 

F=gradG, > (4) 

v=Q v ) 

Bernoulli’s equation for plane unsteady flow can be written 
(from eqs. (14.04) and (9.06) of ref. 4) 

a.+I(u>+e)+^ I -^+lD» ® 

Here a is the speed of sound, and Oq its value in the free 

* In tbo previous literature, both American and British, this is often designated as the 
angle of attack a^sinoe the two are equivalent for a fixed pivot. However, it is 9 but not a in 
•the genera! case of a moving pivot, according to the usual notation of stability analysis. 


stream, where the flow velocity is U. Differentiating this 
expression with respect to time t, and using the fact that 
d[a?l(y— 1)] =a?dplp (ref. 4, eqs. (9.03) and (9.06)) gives 

a t 

— p,=0 (6) 

p 

This, together with the corresponding results obtained by 
differentiating withjrespect to x and y, can be used to elimi- 
nate derivatives of the density from the continuity equation 
(ref. 4, eq. (7.08.2)) ■ 

P/+(p'tt)x+(pP)y— 0 (7) 

The result is that the velocity potential satisfies the equation 

(a 1 — Gx^G^-f- (a 2 — G, 1 ) G„ — 2 G*G/2 — 2 G^G* , — 2 G^Gj, , — 

G„=0 (8a) 

where, from equations (4) and (6), 

a? =(h *+lzl (CP — G* 2 — G„ 2 — 2G ( ) (8b) 

Now introduce a perturbation potential <f>, normalized 
through division by the free-stream velocity U, by setting 

G(r, y, t) = U[x+$(x, y, i)] (9a) 

so that the velocity components are given by 



(9b) 


Then substituting into equations (8) gives 

(1 — M 2 ) $«+$„-2 ^ $„=M 2 [(Y 1 ) 

($x*+$„)+(2<N+4V) 4-xx+*,%,+2 (1+hx) <Iy^+ 

§ ( 10 ) 

For purposes of a second- order solution (and to higher order 
a potential does not exist), the triple products on the right- 
hand side can be disregarded. Thus, the perturbation equa- 
tion becomes finally 

4> w — 0*$xx — 2 ^ 4>«=AP [(7-1) (<!’«+ 

^+ 2 *^+ 2 *^+^ &$«+$,&<)] ( 11 ) 

where — 1. 

PRESSURE RELATION 

Dividing the Bernoulli equation (eq. (5)) by a 2 / (7-1) gives 

£= 1+ ti *.(!-*#«) 02) 
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The flow is isentropic to second order, so that 


y 



and it follows that the pressure coefficient at any point is 
given by 




P~Po_ 

\ptr 


2 

yAP 
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7-1 



(13) 


Substituting the velocity components of equation (9b) , 
expanding in series, and retaining only squares and products 
of perturbation quantities gives finally 


G P =- 2$,— 2 |+a ! -^ J +2W, U +M1 Qff (14) 


Here the second-order solution is required only for evaluating 
the first two terms; the others axe given correct to second 
order by linear theory. 


SMOOTHING OP PROBLEM 

From the leading edge and from any subsequent comer 
of the airfoil spring shock waves or Prandtl-Meyer expansion 
fans that oscillate as the airfoil oscillates. These introduce 
serious complications into the second-order analysis. How- 
ever, the complications can all be circumvented by solving 
a “smoothed” problem in place of the actual problem. The 
solution can thereupon be applied to the actual problem, 
for which it yields the correct result everywhere except near 
the shock waves and Prandtl-Meyer fans. 

The nature of the difficulties can be understood by con- 
sidering first the special problem of steady flow past a single 
wedge (fig. 2). 



The presence of the bow shock wave means that the analysis 
must be undesirably complicated by including the Rankine- 
Hugoniot relations (in a simplified form). A second com- 
plication arises in the differential equation which, for steady 
flow, becomes 

* w - £ 2 <£ i *=M s [(y+1)'M^+ (y- 1)*^ w +2#^] (15) 

In the iteration procedure to be employed, the nonlinear 


right-hand side is evaluated in terms of the first-order 
solution, and the resulting nonhomogeneous wave equation 
solved for the second-order potential. However, for the 
wedge the right-hand side vanishes (to any order), which 
would imply incorrectly that the second-order solution does 
not involve the adiabatic exponent y. More precisely, the 
right-hand side vanishes everywhere except along the Mach 
lines springing from the apex, where it has the singular 
behavior of the Dirac delta function, and only by taking 
account of these troublesome singularities could the correct 
solution be found. 

Both these complications are avoided by the simple device 
of solving the problem of flow past a smooth cusp-nosed 
airfoil of arbitrary shape and then applying the final solu- 
tion to the wedge. It may be imagined that the wedge has 
been smoothed by adding a cusped extension to its nose, as 
indicated in figure 3. It is clear that this artifice removes 



the troublesome singularities from the right-hand side of the 
differential equation. Likewise, it eliminates the need for 
the shock-wave relations because, as indicated in figure 3, 
with sufficient smoothing, shock waves will form only at such 
great distance that their effects cannot reach the airfoil 
surface. Although shock waves axe thus apparently ex- 
cluded, the correct second-order result for the wedge is 
nevertheless recovered from the solution by imagining the 
extension to shrink in size and disappear. The reason is 
that to second order a shock wave is equivalent to the limit 
of a rapid continuous isentropic compression. This limiting 
procedure, which is equivalent simply to applying the solu- 
tion for an arbitrary smooth shape to one that is not smooth, 
yields the proper result except in the vicinity of the shock 
wave (see ref. 5). For an airfoil of general shape, similar 
• broad smoothing must be imagined at any concave corner; 
whereas at convex comers (since no shock waves form) 
the slightest rounding is enough. This smoothing technique 
was applied in reference 4 (p. 399) to steady first-order 
flow past bodies of revolution, and in reference 5 to steady 
second-order plane flow. 

We turn now to the question of generalizing this smoothing 
scheme to an oscillating airfoil. Modification is necessary 
only at the leading edge. Consider first the special case of 
rotation about the leading edge. Then it is enough to con- 
ceive of an extension which is flexible, so that its cusped tip 
can be maintained fixed and directed always into the free 
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stream while the airfoil oscillates, as indicated in figure 
4. (The exact motion of the flexible tip is immaterial, 
provided the surface is sufficiently smooth and its slope re- 
mains small.) After the solution has been found, the flexible 
extension is again imagined to shrink away, and the correct 
result is recovered for the actual airfoil oscillating about its 
nose. 



Figure 4. — Smoothing for airfoil oscillating about leading edge. 


Finally, consider rotation about an arbitrary point. The 
flexible extension must now oscillate in such a way that its 
tip is always directed into the relative wind. Hence, as 
indicated in figure 5, the tip must lie parallel with the free 
stream at the top (and bottom) of each stroke but incline 
in the direction of motion for intermediate positions. ' 


y 



Figure 5. — Smoothing for airfoil oscillating about arbitrary point. 


(i. e., its time rate of change for an observer moving with the 
fluid) vanishes at the surface (see ref. 6), so that 

jSj+H-grad S=0 at #=0 (17) 

With velocity components given in terms of $ by equation 
(9b), and for the smoothed surface described by equation 
(16a), this tangency condition becomes 

$tr=(l+$*)(ef— 0o^ ut gf)~ jjdoe? u, g at y= e f—0 B e {u ‘ g (18) 

where /=/(£), etc. It is convenient to refer this condition 
to the axis y= 0 by expanding in Taylor series. Keeping only 
terms of second order gives 

$j=(l+ $*)(ef— 0oe*“‘ g')—^-0oe ia ‘ g—(ef~ Orff"* g)$n at y = 0 

(19) 

(Here 4> x and on the right-hand side can be evaluated 
from linearized theory.) 

The upstream boundary condition requires that in the 
actual problem, the Rankin-Hugoniot relations (or at least a 
simplified second-order form thereof) be satisfied across an 
oscilliating bow shock wave whoso position must be deter- 
mined. However, shock waves have been eliminated from 
the smoothed problem, so that it is only necessary to require 
that the perturbation potential 4> vanish along the oscillating 
characteristic line (Mach line) springing from the leading 
edge. This insures that all disturbances produced by the. 
airfoil are swept downstream. An equivalent and still simpler 
requirement is that $ and its streamwise derivative vanish 
on, say, the plane 3=0: 


We are accordingly led to consider the motion of an arbi- 
trary flexible oscillating surface described by 

y=ef(x)—0 o e lut g(x) (16a) 

where for the smoothed problem the functions /(r) and g(z) 
have continuous first derivatives. The smoothed problem 
will ultimately be replaced by the actual problem. Accord- 
ing to comparison with equation (3), this means that the 
function g(x) will eventually be identified with ( x—b ). The 
requirement that the leading edge of the smoothed shape be 
always parallel to the relative wind may be written as 


$=$ x =0 at 3=0 (20) 

TRANSFORMATION of perturbation equation 

It is convenient to separate the time-dopendent part of the 
problem from the mean steady flow at zero angle of attack 
(for which the second-order solution is known). Further- 
more, for harmonic oscillations the number of independent 
variables is then reduced to two by separating an exponential 
time factor. Finally, the linear portion of the time-do- 
pendent equation is reduced to normal form by a transforma- 
tion of dependent variable. These threo transformations 
amount to setting 


/( 0)=0 ) 

> (166) 
g'(0)=iwl/U) 

as is clear from equation (19) of the next section. (These 
last conditions, as well as the requirement that j'(x) be 
continuous, must be relaxed in recovering the solution of the 
actual problem.) 

BOUNDARY CONDITIONS 

The boundary condition at the surface of the airfoil is that 

the normal component of velocity is zero. For any surface 

- * 

described by S(x,y,t)=0 moving through a velocity field V, 
this condition means that the substantial derivative of S 


y, t)=4>{x, y)-H 0 e <( “‘- rI) \K3, y) 

where 


__M i « 

K ~Tu 


(21a) 

(21b) 


Here 4> corresponds to the mean steady flow, and the term in 
4 r represents the additional flow associated with the oscilla- 
tion through small angle of <attack. 

Introducing this transformation into the perturbation 
equation (eq. (11)) gives for the potential 4> of the mean 
steady flow 


4>n — P 2 4>xx — Af 2 [ ( 7 — 


( 22 ) 
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whore A is the Laplacian operation 2) 2 /2te 2 +c) 2 /dj/ a ) and for the 
time-dependent part \f r 

<f> e %) I ]-iKA^(y-l)(24^ x +^A^+ 

2 (t+ « t .+fe*+E5 ♦.*.)]- 

^P+Cy— 1 )M*\<j> x gr (23) 

The tangency condition of equation (19) likewise separates 
into the two conditions 

<j> t =e(l+<t> x )f— e<f> t¥ f at y=Q (24) 

^^m/f+ e<tz 4>p,g at y—0 (25) 

For the actual problem the second of these becomes, identi- 
fying g(x) with (x—b), 

+<f, T ) -u j»(x-b)+'(?,-v&)f'- 

£%J+e <tx <l>m(x—b) at y — 0 (26) 

For pressures at the surface of the airfoil, the relation of 
equation (14) can bo expressed in terms of values at y = 0 by 
Taylor series expansion, with the result that to second order 
in thickness and first order in angle of attack 

0 Vu =(,~ 2<fc— 2e<l>xvf+P' i 4>x*—<t>v*)+ 

Moe" *“-«> |jg* *-¥,+« (~ ¥,-¥*)/+ 

P ^x-^v+e^g] (27) 

where all terms are to be evaluated at y=0. 

SOLUTION BY ITERATION 

Although the equation for <j> is nonlinear, that for is 
linear, but with nonconstant coefficients depending upon <t>. 
This corresponds to the physical concept that because of the 
restriction to linear terms in angle of pitch the oscillatory 
part of the flow is an acoustic field with, however, the speed 
of sound varying from point to point in accordance with the 
mean steady flow. 

The well-known linearized or first-order theory results from 
disregarding the right-hand sides of equations (22) and (23). 
Thus, with the first-order potentials denoted by. the lower 
case letters <p and i p, the perturbation equations become 


<Ppp— fP<Pxz=0 

(28) 


(29) 


The second-order solution is obtained by iterating upon 
the first-order results. Using the linear equations to simplify 
the right-hand sides gives for the second-order iteration 
equations 
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4>»- 0‘4> ix =2M 3 [0 2 (N- 1) <p*+<p,% 

and 

l')<px4 / x+ P 2 N<p zx 4'+ fy'Pp] ~ 2A7 <P*4' 

where 


N- 


7+1 M 2 

2 p 


(30) 


(31) 

(32) 


(Here, following the usual subscript notation for derivatives, 
{<p y \ means b<p v 2 /bx, etc.) The second-order solution for 
<£, which leads to Busemann’s well-known result at the air- 
foil surface (ref. 1), was given in reference 5. It is there- 
fore necessary to consider only the second-order problem for 
gf. Details of the interation procedure and discussion of 
its limitations are given for the steady flow in reference 5 and 
apply also to the present problem. 


PARTIAL PARTICULAR INTEGRAL 

The solution of the differential equation for plane or 
axially symmetric steady flow in reference 5 was simplified 
by discovery of a particular integral of the iteration equation 
in terms of the first-order solution. It was also shown there 
that for steady three-dimensional flow a particular integral 
can be found to account for all terms in the iteration 
equation except those involving the adiadatic exponent y 
in the form of N. Likewise, here, a partial particular integral 
that accounts for all terms on the right-hand side of equation 
(31) except those involving N is given by 

S?*=AP(p^) x — (33) 

The complete solution is this partial particular integral plus 
a solution of the reduced equation whose right-hand side 
contains only the terms still unaccounted for: 

gr=2fi 2 N 

W (34) 


FIRST-ORDER SOLUTION 


The first-order solution for <p is known from Ackeret's 
theory to be 

<P=—^f(x—Py) (35) 

It is to be understood here and in all similar expressions to 
follow that this is the potential only for x>fiy, and that <p 
vanishes identically ahead of the bow Mach wave (where 
x<Py). 

The first-order equation for ip (eq. (29)) is most readily 
solved by applying the Laplace transformation with respect 
to x. We denote the Laplace transform of a function either 
by a bar, or by the symbol U, whichever is more convenient 
(and the inverse transform by U -1 ), so that, for example 

UWx))} = Jo e -‘ Zw ® dx (36) 
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Applying this transformation to equation (29), using the fact 
that and vanish, at x=0, gives 

(* 2 +^) ^=0 (37) 

The solution of this equation that represents waves moving 
downstream is 

£=<?(«) (38) 


The coefficient C(s) is determined by the first-order form of 
the tangency condition (eq. (25)), which, transforms to 


where 


?y=— w(s) at y = 0 


w(x)=e ia [g'(x)+vc ^(a;)] 
is the downwash velocity at y= 0. Consequently, 


(39a) 

(39b) 


(40) 


The inverse transformation is readily carried out using the 
standard tables (e. g., ref. 7) together with the convolution 
theorem, which gives as the solution of the smoothed problem 


'/'(V7)=^ "Vo QjV?— /3V) w(x—Qd£ (41) 

The solution for the actual problem is now obtained by setting 
g(x)—x—b, which gives finally 3 

eU *~ e [i +™£i{x-i~b)Jdt 

(42) 

in agreement with the known result of linearized theory (see, 
e. g., ref. 8). 


(eq. (29)), and can therefore be disregarded. Substituting 
the first-order solutions into the right-hand side of equa- 
tion (34) and applying the Laplace transformation gives, to 
order k, 

- ,3 -2eNe-?" [(M s S (x) f (x) } - 

iKS^{f(x)>p(x)}] (44) 

It is readily found that a particular integral of this equation 
is given by an appropriate multiple of ye~^ v . Then, adding 
a complementary function representing downgoing waves 
gives 

*=<7(s) e -^+^y e -^(M’-^ jCU'r } (45) 

where the coefficient 0(s) is to be evaluated from the tan- 
gency condition. Inverting the Laplace transformation 
shows that for the actual problem, in which g(x)=x—b, the 
solution has the form 


,7V 

*=F(z) +^Ly[M*f'+ix(2pzf'- f?bf-f )] (46a) 


where /=/(z). The arbitrary function F{z) is determined 
from the tangency condition of equation (26) to be 


m-Z- 


(2M } — 1 /3 s A , r 

V 2 M* z M s V +e J 


IW-jF-1 


/+ 


p 1 j8"V 2M* " M 2 V 1 'L P 1 

(x-6)/']+i«{[2-(fV- 1) ^~]JV©^+ 

zj+<z-b)f+ 

[^_ (A r- i )] 6/ J (46b) 


SUBFACE PRESSURE COEFFICIENT 


SECOND-ORDER SOLUTION FOR LOW FREQUENCIES 

Because the first-order solution for arbitrary frequencies is 
rather complicated, use is sometimes found for an expansion 
in powers of frequency, which involves only elementary 
functions. The corresponding second-order solution will 
now be carried out in detail, including linear terms in fre- 
quency. This result will serve, for example, to evaluate the 
effects of thickness upon one-degree-of-freedom torsional 
instability, which is primarily a low-frequency effect. 
Thereafter, the result of extending the solution to include 
third powers of frequency will simply be stated. 

POTENTIAL INCLUDING LINEAR TERMS IN FREQUENCY 

Expanding the first-order solution of equation (41) in 
powers of the frequency parameter k and retaining only linear 
terms gives 

^.2/)=^ w(£)d£+0(F)=Tp(z)+ . . . (43) 

where z=x— py. To this order the partial particular integral 
of equation (33) is a solution of the homogeneous equation 


The pressure coefficient at the upper surface of the actual 
airfoil is found from equation (27) to be 




(M*N- 2) 6/']| +O(*V0,0*) 


(47) 


where /=/(x). Here O n is the value for the mean steady 
flow (at zero angle of pitch), which is given by Buscmann’s 
second-order theory. A more useful form of the result is 
obtained by extracting the real part and expressing the 
result in terms of the instantaneous angle of pitch 0(f) and 
its time rate of change 0(f). Furthermore, the parameter e 
has served its purpose of distinguishing terms of different 
orders and can be eliminated (according to eq. (1)) . Thus, 
on the upper surface of an arbitrary airfoil that is described 
at zero angle of pitch by y=Y(x), is pivoted about a point a 
distance b downstream of its leading edge, and performs 
slow angular oscillations described by 0(f), the pressure 


* The smoothing was, ol course, unnecessary in this first-order problem. 



SUPERSONIC FLOW PAST OSCILLATING AIRFOILS INCLUDING NONLINEAR THICKNESS EFFECTS 


617 


coefficient is, to second order in thickness and first order in 
angle of pitch, 


0,=0„-je +: 

2 p W- i) F . 


(2 -M*) (MW— 1) 


xY'+ 


MW- 2 
ft 


6F']^+O(0, e s 0,ft) 


(48a) 


u 



Figure 6. — Velocity relative to leading edge of oscillating airfoil. 


Here the value for the mean steady flow is (ref. 1) 



MW-2 

ft 


Y /2 


(48b) 


(A preliminary report of this result was given in ref. 9.) In 
this form, the result is not restricted to sinusoidal motion 
but applies to any oscillation that is sufficiently smooth and 
slow that the pressures depend significantly only upon the 
instantaneous angle of pitch and angular velocity. 

The pressure on the lower surface of the airfoil is obtained 
from these equations by reversing the sign of 8, and taking 
F( x) to be the ordinate of the lower surface, measured 
positive downward. 

The result for plunging motions can be extracted by 
letting 6 tend to zero and b tend to infinity in such a way 
that their product remains finite, say 

bd(t) = -h(t) = -l k> e i “ t (49) 


In the limit, the airfoil simply translates vertically according 
to y——h(t). The pressure coefficient on the upper surface 


is 


4 


C,=C,-- a T - 2 


2 n „ MW — 2 


P u 


ft 


y, _ 

1 u 


(50) 


CHECKS ON THE RESULT 

The solution can be tested in several special cases for 
which the result can be derived from other considerations. 

Of the five terms in equation (48a), the first is known 
from Busemann’s steady second-order solution, and the 
second and third from linearized unsteady theory. The 
fourth is obtained by using the instantaneous airfoil slope 
(Y'—d) instead of the mean steady slope Y' in Busemann’s 
formula and retaining only linear terms in 9. Therefore, 
only the last term, which is the essentially new result of the 
present analysis, requires verification. 

Just at the nose of an oscillating airfoil, the pressure can 
be determined exactly if the bow shock wave is attached. 
The transition through the moving bow shock is instan- 
taneous, and so depends only upon the relative velocity at 
that instant (see ref. 4, p. 297). Hence the pressure just at 
the nose is instantaneously the same as on a wedge of the 
same vertex angle in steady flow with the same relative 
velocity. In the present problem, the relative velocity is 
compounded of the horizontal velocity U of the free stream 


< Note that, as It should be, this is Just the result ol using Bnsemann’B formula for steady 
Dow (eq. (48b)), with the local slope decreased by the Instantaneous apparent angle 1/17; 
see the discussion In the following section of conditions at the nose. 


plus the instantaneous vertical velocity of the leading edge, 
which is given by 6b (see fig. 6). The effect of the vertical 
component upon the equivalent free-stream velocity and 
Mach number is of second order in angle of pitch, but the 
equivalent vertex angle of the airfoil is increased by the 
apparent downwash angle 6b/U. Replacing Y' by F'+ 
6b /U in Busemann’s formula (eq. (48b)) gives, to first 
order in angle of pitch 


°'.-°:+l i T} +2¥ ^W AlT 'V (51) 


which checks the part proportional to b of the last term in 
equation (48a). 

The remainder of the term in question can be checked for a 
single-wedge airfoil oscillating about its vertex (fig. 7). It 


u,M 

Po’Po 



Figure 7 . — Wedge oscillating about its vertex. 


can be shown using the results of reference 3 that in this 
case disturbances reflected-from the shock wave are of third 
order in the wedge angle (although for other pivot positions 
they are of second order). Therefore, a solution correct to 
second order in thickness and first order in angle of pitch 
can be found by applying linearized theory to the mean 
steady flow behind the shock wave.® For slow oscillations, 
the first three terms of equation (48a) give 


V-V i 

i 


ft 


2(2— Mi 2 ) _ 6 

'3 Xl 


ft 3 


Ui 


(52) 


* This concept was suggested to the author by W. P. 1 ones of the National Physical Labo- 
ratory, England. 
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where subscript 1 denotes values in the mean steady flow 
behind the shock wave. From linearized theoiy 


Mi=M [1 — P(N — l)e] ' 

Pi=P [l— " (iV— l)ej 

u '- u H) 

Pi=Po (l+Vjf- e) 

<"=»( 1 +T < ) 


(53) 


where e is the semivertex angle. Hence, referring the sur- 
face pressure coefficient to free-stream quantities (and noting 
that Zi—x to second order) gives 


0*-0, o-| e+2 ljj£ z 1-2 <0+ 


2 2M\N- l)+(2-M 2 )(MW- 1) ^ ^ m 


which checks the first two parts of the last term in equation 
(48a) when Y (z) — ex. 

Recently, Idghthill has given a further check for the case 
of Mach numbers so high that 1/M 2 is negligible compared 
with unity (ref. 10). 

LIFT AND MOMENT COEFFICIENTS FOR SYMMETRICAL AIRFOILS 

The coefficients of lift and pitching moment (about the 
pivot) are given in terms of the pressure coefficients on the 
upper and lower surfaces by 


4T/2-M* 2M 2 — 3 b b*\ M 2 (A 7 -1) b ,, 
pl\~W h 2/3* c c 2 / c * 

jj- Tdx+2 ^=lj£^ b)Tdx -]§ m 

Inside each square bracket the fust term is the result of 
linearized theory, and the remainder represents the second- 
order effect of thickness. Thickness effects are seen to 
appear in the form of the area of the airfoil profile and its 
first moment about the vertical line through the pivot. 

If the airfoil has a blunt trailing edge of semithickness 
Y(c), the following additional second-order terms must be 
added to the above expressions: 
to Ci'. 

4 MW— 2 Y(c) 4 TMW-2 b (2-M 2 )(MW-l)l Y(c) c6 

P P c p [_ P c + P 3 J c V 

(69) 

to c m : 

(60) 

EXAMPLES: BICONVEX AND DOUBLE- WEDGE AIRFOILS 

To second order, a biconvex airfoil of thickness ratio r is 
given by the parabolas 

V=±T(x), Y(x)=2r^(c-x) (61) 


(C P —C Pu )dx (55a) 

Cm= (?Jo ib—z)(C P —C Pu ) dx (55b) 

For simplicity, consider only airfoils symmetric about the 
chord line. In this case, the pressure difference is given by 


MW- 2 


0,-0, x+b^l+i —p 


Y'e- 


4 pM 2 W-l) Y ] (2 —AI~){M 2 N— 1) xT , + 


MW — 2 


bY' 


U 


(56) 


If the airfoil has a sharp trailing edge, substituting into 
equations (55) and integrating by parts gives 




'2 — M 2 , b\ 

, MW-3M s +2 1 

. 2/S 2 1 c) 

1 p 3 c 2 


J>*]u < 67 > 


and 


The expressions for lift and moment become 
/ '2—M i , b 




r/ 2 -M 2 , b\ 

A 2/S 2 1 c) 
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2/3 3 


MW — 2 
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(62) 


5) 


+ 


(63) 


It happens that for a double-wedge airfoil (with maximum 
thickness at midchord), both the area and first moment are 
just three fourths of those for the biconvex airfoil. Conse- 
quently, the above results apply to double-wedge airfoils if 
t/ 3 is replaced throughout by r/4. 

In the expression for pitching moment, the term propor- 
tional to 0 represents an aerodynamic stiffness or restoring 
moment in phase with the angular displacement, while the 
term proportional to 6 corresponds to an aerodynamic damp- 
ing moment in phase with the angular velocity. The effects 
of thickness upon aerodynamic restoring and damping 
moments are shown in figures 8 and 9 for a 6-percent-thick 
double-wedge or a 4}£-percent-thick biconvex airfoil with 
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Fioure 8. — Effect of thickness upon restoring-moment coefficient for 
4)4-percent-thick biconvex or 6-percent-thick double-wedge airfoiL 



Figure 9. — Effect of thickness upon damping-moment coefficient for 
41 ^-percent-tbick biconvex or 6-percent-thick double-wedge airfoiL 


two different pivot positions. The figures have been labeled 
with both the usual American notation (e. g., ref. 11) and the 
British notation (e. g., ref. 12) for flutter derivatives, defined 
for harmonic oscillations of arbitrary frequency by 

— i (American) 

A 

~ (m e +i\me) (British) 

z 

Here k and X are the reduced frequency in the American and 
British notations, respectively, related by 




X=2k=% 


( 65 ) 


For slow oscillations, the flutter coefficients are obtained 
from equation (63) according to 




mi—- 


2 Mi ~2 c bd ' 


( 66 ) 


(and, indeed, this was originally the definition of me and mi, 
that given by equation (64) being a later extension to the 
case of rapid harmonic oscillations) . 

It should be noted that according to second-order theory 
the nonlinear effects of thickness are themselves linear in 
thickness. This means, for example, that doubling the air- 
foil thickness ratio would double the distance between the 
linearized and second-order curves of figures 8 and 9. 


NEUTRAL DAMPING BOUNDARY 

Linearized theory indicates the possibility of instability 
of pitching oscillations for low frequencies. For a range of 
Mach numbers below ^/5/2= 1.58 and pivots ahead of two- 


thirds of the chord, the aerodynamic damping moment 
becomes negative, and so tends to destabilize. (Whether 
or not the motion is actually unstable depends, of course, 
upon the other dynamic parameters in the problem.) This 
zone of possible instability shrinks and eventually disappears 
as the frequency of oscillation increases. The present low- 
frequency solution is therefore adequate for determining 
how the region of instability is modified by nonlinear thick- 
ness effects. 

Figure 10 shows the boundary of neutral aerodynamic 
damping for slow oscillations of a 4 %-p ercent-thick biconvex 
or 6-percent-thick double-wedge airfoil. The aerodynamic 
damping is destabilizing for Mach numbers and pivot posi- 
tions lying inside the loops. Within the region where linear- 
ized theory predicts a destabilizing moment, thickness is 
seen to exert a further destabilizing effect except for pivots 
near midchord. The second-order solution becomes unreli- 
able when the bow shock wave detaches, at about M= 1.2. 

COMPARISON WITH PREVIOUS INVESTIGATIONS 

Two previous investigators have sought a second-order 
solution for slowly oscillating airfoils in supersonic flow. 
Their results agree neither with each other nor with the 
present solution. 

In 1947, W. P. Jones obtained an estimate of the thickness 
effect by assuming that the ratio of second-order to linearized 
pressure disturbances is the same for slow oscillations as 
that given by Busemann’s formula for steady flow (ref. 13). 
That this assumption is not altogether correct is indicated 
by the fact that the results do not check those obtained for 
a wedge oscillating about its vertex by applying linearized 
theory to the mean steady flow behind the shock wave. 
However, the assumption is correct at the leading edge, and 
also (as noted by Ldghthill in ref. 10) in the limit of high 
Mach number. It is seen in figure 10 that this estimate 
fails to give a useful prediction of the actual effects of thick- 
ness, except for pivots near midchord. 

In 1951, Alexander Wylly attacked the problem by 
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Figure 10. — Neutral damping boundary for 434-percent-thick biconvex 
or 6-percent-thick double-wedge airfoiL 
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methods similar to those used here (ref. 14). Unfortunately, 
it appears that the smoothing was not carried out with 
sufficient care; as a consequence, the solution satisfies none 
of the three checks discussed previously. In contrast to 
the present results, the effect of thickness upon aerodynamic 
damping was predicted to be stabilizing and so great that 
for airfoils of the thicknesses shown in figure 10 the zone of 
possible instability would have disappeared altogether. 

Concurrently with the original appearance of the present 
work, Martin and Gerber have published an independent 
investigation of the second-order effects of thickness on the 
stability derivatives for airfoils in constant vertical accelera- 
tion and constant pitching (refs. 15 and 16). Their results 
agree completely with those given above. 

COMPARISON WITH EXPERIMENT 

The first- and second-order theories are compared in 
figure 11 with the results of recent English experiments car- 
ried out at the National Physical Laboratory. 8 The aero-, 
dynamic damping at low frequencies was measured for 
biconvex airfoils 5 and 7.5 percent thick at Mach numbers of 
1.42 and 1.61 with various pivot positions. For the points 
shown the amplitude of angular oscillation was 1.5°; increas- 
ing the amplitude to 3° was found to affect the results only 
at the lower Mach number. Including second-order thickness 



Pivot position, b/c 
(a) 711=1.42 


Figure 11. — Comparison of experimental and theoretical damping- 
moment. coefficients for biconvex airfoils. 

• The author is Indebted to W. P. Jones, J. B. Bratt, and W. E. A. Acorn of the National 
Physical Laboratory for kindly making these data available in advance of publication. 
The work is to bo published under the title “Measurements of Pitching Moment Dertva- 
Uves for Two-Dimensional Models at Subsonio and Supersonic Speeds, and for a Rectan- 
gular Model of Aspect Ratio I at Subsonio Speeds,” by Bratt, Raymer, and Townsend. 
(It may be noted that earlier experimental results from the National Physical Laboratory 
that were cited in a previous version of this paper (NACA TN 2932) are now believed by the 
oxperimentets to be unreliable.) - ' 



(b) M=1.61 
Figure II. — Concluded. 


effects in the theory is seen generally to improve the agree- 
ment with experiment. The lower Mach number is close to 
the limit of purely supersonic flow — M= 1.38 for tlio 7.6- 
percent-thick section — at which the theory presumably 
breaks down. 

EXTENSION TO CUBE OF FREQUENCY 

The dependence of thickness effects upon frequency of 
oscillation can be estimated by extending the second-order 
solution to include higher powers of frequency. This has 
been carried out for an arbitrary airfoil by including second 
and third powers, which is enough to show an effect of 
frequency upon both aerodynamic stiffness and damping. 
The computation, though cumbersome, is a straightforward 
extension of the previous analysis, so that only the final 
result will be given here. 

The expression for pressure coefficient on the upper surface, 
corresponding to equation (47), is found to bo 
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. , T(17M 4 - mP-4:)N-(5M 2 -2)(4M 2 - 1) , 

2M 4 /3* J 

UP-(5AP-2)N , f, , 

2M 4 J 

7M\2AP- 1) - (12M 4 -3M 2 -4)AI , 

2M 4 /3 2 X J 

2(AP+l)-(AP+S)N ,, , 


(3M 2 +2)-(3M 2 +4)A r ^ , 
4M 4 /3 2 


2(M 2 + 1) — 3MW 




5M 2 +2-M 2 (4+M 2 )AT 

12W ^ J (6?) 

where/ s/(a) and J/ = J" /(f) d£, etc. (The first four terms are 

the result of linearized theory.) This result meets the three 
tests discussed previously, and also checks the solution given 
later for a wedge with general pivot location at arbitrary 
frequency. The resulting expressions for lift and moment 
involve the airfoil thickness in the form of the area of the 
profile and its first three moments about the vertical line 
through the pivot. 

EXAMPLE: BICONVEX AIRFOIL 

These rather formidable results simplify considerably for 
specific airfoils. For example, for a biconvex airfoil of thick- 
ness ratio r oscillating about midchord, the pitching-moment 
coefficient is given by 

O m 2—AP 2 AP APm*-l) , MW— 2 

40oe iut 12/3 3 ■** 16/3 5 160^ 7 “ t " T 3/S 2 

MW- 1) , 12(2— M 2 )+(M 2 — 4)iV 

tXr +X tM 2400 s + 


iWAP 


2(48M 4 — 61 M 2 +l)-(37M 4 +25M 2 -16)iV 
1440/S 8 


( 68 ) 


The first three terms are the result of linearized theory. 

The component of this moment that is out of phase with 
the angle of pitch gives a parabolic approximation for the 
variation of aerodynamic damping with frequency as shown 
by dashed lines in figure 12 for a 5-percent-thick biconvex 
airfoil. The accuracy of the parabolic approximation for 
linearized theory is indicated by comparing it with the exact 
result (solid line) . In this example the linearized and second- 
order curves run almost parallel, which means that the non- 
linear effects of thickness vary only slightly with frequency. 

Recently, Jones and Skan have treated biconvex airfoils 
at arbitrary frequency by a numerical procedure (ref. 17). 
Their result is shown in figure 12. It fails to give the initially 
parabolic form of the curve that is implied by the fact that 
the second-order solution, like the first-order result, can be 
expanded in powers of the square of the frequency. Their 



Reduced frequency, \*Zk -ac/U 

Figure 12. — Effect of frequency upon damping-moment coefficient for 
5-percent-thick biconvex airfoil pivoted at midchord. 


solution involves several doubtful assumptions, in particular, 
that the effect of the bow shock wave can be disregarded. 
It has already been remarked here that, actually, the bow 
shock has a second-order effect unless the pivot lies at the 
leading edge. 

STABILITY DERIVATIVES 

The oscillating motion of the airfoil has heretofore been 
described in terms of its angle of pitch 0 and (negative) 
elevation h, as is customary in flutter analysis, and the 
flutter derivatives have been calculated. In stability analy- 
sis an alternative pair of coordinates is usually employed: 
the angle of attack a with respect to the relative wind and 
the rate of pitching 2=0. The motion is then no longer 
restricted to harmonic oscillations, so that there are actually 
an infinite number of stability derivatives. However, only 
the first three — c„ a , c mq , — for moment, and their counter- 

parts for lift, are ordinarily considered significant. 

Steady-flow theory gives c„ a . It is shown in Appendix 
B (eq. (B12)) that the combination (c mj +c„ A ) is given by the 
solution for low frequencies, but that c„ Q and c m& separately 
can be found only from the solution that includes the square 
of the frequency. For this purpose, only the solution for 
plunging is required, and this can be extracted from equation 
(67) as before by letting the pivot recede to infinity and the 
angular amplitude diminish according to equation (49), 
which gives 

C,~C Pl =^ e P^/^+W-i) xf (z)]J 

(69) 

Integrating for the moment according to equation (55b), 
replacing/ by 7 according to equation (1), and then extract- 
ing c m& according to equation (B12) gives 

The combination (c m? -f- c„ 6 ) is given, according to equation 




622 


REPORT 1183 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


(B12), by twice the coefficient of (c6/U) in equations (58) 
and (60). Together with equation (70), this gives 

<?i> 

These results have been derived independently by Martin 
and Gerber (refs. 15 and 16). The agreement serves as a 
further partial check on equation (67). 

SECOND-ORDER SOLUTION FOR ARBITRARY FREQUENCIES 

For some purposes the previous solution for slow oscilla- 
tions may be inadequate. In principle, the second-order 
solution can be extended to include still higher powers of 
frequency, but the labor required is clearly prohibitive. 
Alternatively, one can attack directly the problem for arbi- 
trary frequencies. 

The second-order solution can, in fact, be carried out for a 
general airfoil at arbitrary frequencies. However, the result 
is formidable, involving multiple integrals of products of 
Bessel functions, and the reduction to simpler form for 
specific profiles appears to be difficult. 

A more practical approach is to choose a specific airfoil 
shape in advance. Then the second-order solution involves 
only functions of the type encountered in the linearized theory. 
In particular, it is found that (at least for the simplest shape) 
the final expressions for lift and moment involve only func- 
tions which have been already studied and tabulated, so that 
numerical results are readily obtained. 

MODIFIED SMOOTHING PEOCEDUEE 

The smoothing discussed previously must be dropped at 
an earlier stage of the solution when a specific airfoil shape 
is chosen. It is therefore necessary first to modify the 
differential equation and boundary conditions so that no 
singular terms appear. 

Consider first the differential equation (eq. (31)). Apply- 
ing the Laplace transformation of equation (36) , and envision- 
ing the smoothed problem, so that ip and fa vanish at 2 = 0 , 
reduces it to 

(f+jp) ¥=2 (M*8-uc) [w-1)C{^} + 

-C{ <Ptfa } — w -C{ <PxP}~^ (72) 

In this form, all troublesome second derivatives have dis- 
appeared from the right-hand side, so that the smoothing 
can now be dropped. 

Consider next the tangency condition of equation (25), 
which contains second derivatives of both <p and- p. These 
cause no difficulty when, as in the solution for slow oscilla- 
tions, the tangency condition is imposed as it stands. How- 
ever, its Laplace transform will be used here, and then the 
Dirac delta functions associated with the second derivatives 
would affect the integration implied in the inversion of the 
transformation. These troublesome second derivatives can 


be eliminated by first expressing them in terms of x deriva- 
tives only through the first-order equations : 

<Pm=&‘<Pxx 
'/'w=|5 s 

and then applying the Laplace transformation while envi- 
sioning the smoothed problem (so that <p, ip, fa, and <p x vanish 
at 2 = 0 ), which gives 

%= — io(s) { e itc <p x ] -feAftC { PJ ' } — {'Pf'} — 

/9 s (s— fic)-C{e* I Vx(z— &)} aty=0 (74) 

Again, all second derivatives have disappeared, so that the 
smoo thing can be dropped. 

Consider finally the upstream condition. For the smoothed 
problem, one statement of the condition was seen to be that 
the solution represents downgoing waves. This moans that 
the complementary function for the iteration equation should 
have the same form as the linearized solution, and this state- 
ment of the upstream condition applies as well to the actual 
problem. 

SOLUTION FOE WEDGE 

The simple case of a single wedge illustrates the method of 
solution that can, in principle, be applied to any profile 
formed of piecewise analytic arcs. For a wedge of semivertox 
angle e, the first-order solution for <p is given by equation (36) 
as 

<P=~ (*-/ By) (76) 

and the first-order solution for p is given by equation (42). 
Substituting into the right-hand side of equation (72) gives 
for the transformed iteration equation 7 

0 * 2 «(M s a— fic) |~1+ 

(76) 

A particular integral of this equation is given by an appro- 
priate multiple of and adding a complemen- 

tary function that represents downgoing waves gives tho 
solution 

?= (7(8)e-*V s! +*’Mfi + 

1 rAfto—tA T I , (N—IJs—ixN/M 2 '} 

^ ) l-,tf+7jW + W/M* ] w[s)ye 

(77) 

The coefficient (7(s) is evaluated by imposing the tangency 

* In this simple case, the partial particular integral of equation (33) offers little actual sim- 
plification. 
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condition of equation (74), with the result that the Laplace 
transform of the second-order solution is found to be 


__ c~ Po V « , +* a /-W 3 / r 

i (m-i,) w-w-ww a + 

P Vs’+^/M 3 V 


l+€2/(M*s— tK)+ 


(78) 


The inversion can be carried out using the standard tables 
(e. g., ref. 7) together with the convolution theorem. For 
calculating surface pressures, it suffices to obtain the solution 
in the plane y= 0, which is found to be 


'P(z,0)==ij^ J 0 (jjt) w(x—£)d£+ ex w(x)+ 

in(2N— 1 )w(x— £)~ k* e U(x ~Q 


whore 


w(x)=e ia ^1+iK jp fa— 6)J 


With the pivot at the nose (6=0), this agrees with the result 
of applying linearized theory to the mean steady flow behind 
the shock wave. Also, when expanded in powers of fre- 
quency, it agrees with the previous low-frequency solution' 
up to terms in k 3 . 

The surface pressure coefficient can now be calculated 
from equation (27) and the lift and moment coefficients 
from equations (55). 

EXAMPLE: WEDGE PIVOTED AT NOSE 

For simplicity, the results will be given only for the special 
case of rotation about the nose. Then it is found that the 
lift and moment coefficients are given by 
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where v=0 for the lift and v= 1 for the moment. The func- 
tions/,, given by 

/„(M,X)= (M Xz) dx (81) 

arise in linearized theory for n ranging from 0 to 3. They 
have been studied and tabulated by von Borb61y (ref. 18), 
Schwarz (ref. 19), and Garrick and Rubinow (ref. 11). They 
can all be expressed in terms of / 0 by a recurrence relation 
due to von Borb61y (see, e. g., eq. (A.87) of ref. 12), so that 
the additional / 4 required here is easily computed. 



Reduced frequency, X*2 k*u>c/U 


Figure 13. — Damping-moment coefficient for 10-percent-thick wedge 
oscillating about its vertex at M=10/7. 

Figure 13 shows the variation of aerodynamic damping 
moment, according to first- and second-order theories, for a 
10-percent-thick single wedge oscillating about its nose at a 
Mach number of 10/7. Also shown for comparison are the 
parabolic approximations of the low-frequency analysis. It 
is seen that the thickness effect is reversed at high flutter 
frequencies, as is suggested by the parabolic approximation. 

DISCUSSION 
HlGH£R*OHD£B EFFECTS 

The moderate magnitude of second-order effects would 
suggest that the influence of third- and higher-order terms is 
of no practical importance, except perhaps in the transonic 



Pivot position , b/c 


Figure 14. — Neutral damping boundary for wedge of 5° semi vertex 

angle. 
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range near shock detachment. This supposition can be con- 
firmed in the case of the single-wedge airfoil, for which a 
solution exact in thickness (but linearized with respect to 
angle of pitch) has been derived in references 2 and 3 . Figure 
14 compares the boundaries of neutral aerodynamic damping 
for a slowly oscillating wedge of 5° semivertex angle as pre- 
dicted by the linearized, second-order, and exact theories. 
The second-order solution lies close to the exact result down 
to the Mach number for shock detachment (which is almost 
the same as the Mach number at which the flow ceases to be 
purely supersonic). 8 

APPLICATION TO FINITE-SPAN WINGS 

Extension of the second-order solution to wings of finite 
aspect ratio does not seem possible at present. No second- 
order solution has yet been found even for steady flow past 
the simplest lifting wing. 

Fortunately, the main conclusion to be drawn from the 
present analysis is that nonlinear thickness effects are quite 
moderate in magnitude. Practical supersonic wings will, 
therefore, probably be so thin that nonlinear effects are neg- 
ligible, so that reliance can be placed in the predictions of 
linearized theory. Only if the wing is unduly thick, or if 
the Mach number is close to unity, or if unusual accuracy is 
required, may the engineer be forced to estimate the effects 
of thickness. In this event, he might assume that the effects 
of thickness are in some sense additive to those of aspect 
ratio, provided the aspect ratio is high and the frequency low. - 
For example, the two-dimensional correction might be ap- 
plied stripwise to the spanwise loading predicted by linearized 
theory. Some indication of the extent to which such an as- 
sumption would be justified can be obtained by considering 
other pairs of effects whose combined influence is known. 
Figure 12 shows that the effects of thickness and frequency 
are roughly additive for the frequencies of usual practical 
interest in dynamic stability analysis (say \<0.2). Like- 
wise, figure 15 shows that the effects of aspect ratio and fre- 
quency, determined from Watkins’ linearized solution for the 

• In flguro 14 the curve designated second-order theory includes all second-order effects 
and also some (but not all) of the third-order effects. It therefore does not represent the 
second-order solution of this report. The curve for that solution (calculated from eqs. (68) 
and (60)) can be plotted from the following values: 

6/e: -L6 — LO —0.5 0 0.33 0.6 0.0 

M: 1.21 L26 L34 L52 1.63 1.62 1.26 

The revised curve faTU significantly below that shown for the exact theory. Henoe the 
conclusion reached in the discussion of higher-order effects on pages 623 and 624 most be modi- . 
fled; third- and higher-order terms can be of some practical significance. The author Is 
indebted to Dr. Georg Drougge for pointing out this discrepancy. 
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Figure 15. — Effect of aspect ratio upon damping-moment coefficient 
for rectangular wing according to linearizod theory. 


rectangular wing (ref. 20), axe nearly additive in the same 
range of frequencies. 

Martin and Gerber have calculated the second-order effects 
of thickness upon damping in roll for wings of infinite span 
(ref. 21). They then estimate the effect for a finite rectan- 
gular wing by increasing the result of linearized theory in 
the same ratio as for the infinite wing. Agreement with ex- 
periment is thus considerably improved, which gives further 
assurance that superposition of thickness effects may bo justi- 
fied also for oscillating wings. 

Recently, Acum has estimated the aerodynamic damping 
of pitching oscillations of rectangular wings at supersonic 
speeds by assuming that the thickness effects of the presont 
theory can he added to the aspect-ratio effects predicted by 
linearized theory (ref. 22). 

FURTHER ANALYSIS 

If thickness effects of the magnitude indicated by the 
present analysis are judged to be significant at flutter fre- 
quencies, extension of the high-frequency solution to more 
practical profiles would be warranted. The solution given 
for a single wedge at arbitrary frequency should he extended 
next to the biconvex or double-wedgo airfoil. Although con- 
siderable computation is involved, it does not appear that 
the labor would be prohibitive. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. SO, 1953 




APPENDIX A 

SYMBOLS 


a 

A 

b 


C(s) 

Cl 
C m 

Cm a , C m ) 

Cm & ) 
Cmg> CmU 
c mj) Cm/,) 
0, 

0 v o 

n n 

1/ pi 


m 

fn 

y( x ) 

hit) 

ho 

1c 

m h , m e> ) 
ml, mi ) 
M 

Mi, M a ,\ 

Me, Mi ] 

N 

P 

2 

8 

S(x, y, t) 
t 

u, v 


U 

V 

wix) 


speed of sound 
aspect ratio of wing 

downstream distance from leading edge to pivot 

airfoil chord 

constant of integration 

section lift coefficient 

section moment coefficient 

stability derivatives (See eq. (B5) .) 

flutter derivatives (See eq. (B6).) 
pressure coefficient 

surface pressure coefficient in mean steady flow 
pressure coefficients on upper and lower surfaces 
of airfoil, respectively 

function defining upper surfaco of airfoil at zero 
angle of pitch 
See equation (81). 

function defining amplitude of oscillation of 
upper surface of smooth airfoil 
elevation of plunging airfoil, positivo downward 
amplitude of harmonic plunging oscillation 

reduced frequency, *■ 

British flutter coefficients (See eq. (B9) .) 

free-stream Mach number 

American flutter coefficients (See eq. (B8) .) 

7+1 M 3 
2 j9 3 

static pressure 

rate of pitching, 6 

Laplace transformation variable 

function defining moving surface 

time 

velocity components parallel and perpendicular 
to free stream 
speed of flight 

velocity vector 
downwash velocity at y = 0 


x,V 


Tix) 

z 

a 

|S 

7 

A 

€ 

e 

e 0 


x 

V 

p 

T 

<P 

<t> 

$ 

* 

** 

W 

a 

( )o 
( )i 

(•). 
O ] 
•C{ M 


coordinates parallel and perpendicular to flight 
direction, moving with mean steady velocity 
of airfoil 

' ordinate of upper surface of airfoil at zero angle 
of pitch 
x—i 3y 

angle of attack — angle of airfoil with ’respect to 
relative wind at pivot 

Vm^i 

adiabatic exponent of gas 

Laplaeian operator, j^ 2+^“2 

small parameter representative of airfoil thick- 
ness 

angle of pitch — angle of airfoil with respect to x 
axis 

angular amplitude of harmonic pitching oscilla- 
tion 
M^co 

jS V 

reduced frequency, 

constant that is 0 or 1 in equation (80) 
density 

airfoil thickness ratio 

first-order mean steady perturbation potential 
second-order mean steady perturbation potential 
complete perturbation potential . 
first-order time-dependent perturbation poten- 
tial 

second-order time-dependent perturbation poten- 
tial 

partial particular integral of time-dependent 
iteration equation 
angular frequency of oscillation 
complete velocity potential 
free-stream conditions 

conditions in mean steady flow behind bow shock 
wave on wedge 

differentiation with respect to time 
Laplace transform 
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APPENDIX B 

CONNECTION BETWEEN FLUTTER AND STABILITY DERIVATIVES 


la flutter analysis the motion of the airfoil is described by 1 
the angle of pitch 9 and the (negative) elevation h, whereas 
in stability analysis it is described by the angle of attack a 
and the rate of pitching g. From figure 16 it is seen that 


tan a 


tan d+K/U 
1 — (fi/U) tan 6 


CB1) 



Figtoh 16. — Alternative coordinates for describing motion of oscillating 
airfoil through still air. 

and by definition g=0. Hence to second order in the angles, 
these alternative coordinate systems are related by 

a=6+ ^\ CB2) 

a=* ) 

For a given profile, gas, Mach number, and pivot, the 
moment coefficient at any instant depends on the entire 
previous history of the airfoil motion. Thus in stability 
coordinates 

C m (t)=C m L a (f), am (B3) 

Here the heavy brackets are used to emphasize that this 
is not an ordinary function, but a functional of the two 
coordinates, which are themselves functions of past 
time. Now if a and g are analytic functions (for all past 
time), they can be resolved into their successive time deriva- 
tives. Then c m assumes the nature of an ordinary function of 
these infinitely many variables, each of which depends upon 
time as a parameter: 

c»(f)=c„[o(f), a(t), ■ ■ ■ , 2 (0, 2 (0, • - • ] (B4) 

Finally, if c m is analytic in each of these variables, Taylor 
series expansion gives 3 

i British writers use r instead of h. Both American and British writers have previously 
used a Instead of e t but 6 la preferable for the r ea s on given in footnote 2, and is apparently 
now being adopted by the British. 

a The reader may choose to regard this result as being obvious from physical considerations. 
In which case the preceding mathematical formalism Is superfluous. 
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. . . (B6> 

where only linear terms have been retained. Here the 
coefficients have been made dimensionless using the charac- 
teristic time c/2 U, in accordance with the usual American 
notation. In stability analysis the motion is usually assumed 
so slow and smooth that the three terms shown explicitly 
are sufficient. The coefficients c„ a , c„ t , c m/t> etc., which are 
(aside from factors of c/2U) the first partial derivatives 
of the function c m , are the stability derivatives. 

Proceeding similarly in flutter coordinates would lead to a 
corresponding expansion in terms of flutter derivatives: 

C m ({) = d(f) 6 (t)c md +jj (^ 2 / 7 ) • • • 

(B6> 

Here the first term has been omitted since it is clear physi- 
cally that c„ h vanishes identically — the moment depends 
on changes in elevation but not on the elevation itself. 
However, it is customary in flutter analysis to consider only 
harmonic oscillations: 

d(f)=0 o e*"‘,h(t)=he? ut (B7) 

and then the infinite series for c m collapses to just four 
terms. In the usual American notation for flutter coefficients 
(ref. 11) 

<v=-: 2P (2 \ M,+dM 3 +^j (B8> 

and in the usual British notation (refs. 12, 13), with ( 6 , h ) 
in place of (a, z ) 

c*= 2 0 m h +^ m- h +em e +jj (B9) 

Comparing these expressions with equation (B6) in the 
special case of harmonic oscillations shows that 

— 2 ¥ Mi = 7714= — . . . 

2 7c Af2=277l4=Cn,4 ' * * 

—2Je 2 M 3 =2mg=c ai —Jc 3 c^+ . . . 

— 2k M Al =\mb=c*$— . . . 

The relations between derivatives in the two coordinate 
systems are found by regarding each system as consisting 
of an infinite number of coordinates that are related by 
equation (B2) together with all the relations obtainod by 
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differentiating equation (B2) repeatedly -with, respect to 
time. Hence 

Cm$ — Cm a C*t a =C»ij 

h = Cm a Cm t =C m £ CmJ 

= C ma = C m ]J 

e »h~ C m a 


Finally, combining these relations with those of equation 
(B9) gives the omnibus relations 

— 2Jc 1 2 M 1 = m b =— A 3 4 * 6 c„ x + . 

— 2k Mi=2mh=c m - h ~ ... 

—2k i M a =2me=c mB — . . . 

— 2k M t =4:m=c„j— . . . 

The corresponding relations for lift can be found in the 
same way. 

Equation (B12) shows that the combination (c„ + 
c„ A ) is proportional to the first term in the expansion 
of the flutter coefficient rm or in powers of frequency. 
It is therefore given theoretically by the solution linear in 
frequency, and experimentally by measurements at low 
frequency. However, to determine c„ & alone, and hence 
c M(I , squares of frequency must be retained. For this 
reason, it is difficult to find c„ s and by oscillating an 
airfoil in a conventional wind tunnel. It would be necessary 
to perform the experiment at various frequencies and so to 
determine the initially parabolic variation with frequency 
of the component of moment that is out of phase with the 
angle of pitch. 
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